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SPIN-ORBIT COUPLING AND ENERGY SHIFTS IN SINGLE 


CRYSTAL AND PYROLYTIC GRAPHITE 
by John A. Woollam 
Lewis Research Center 

SUMMARY 

A study of the spin-orbit coupling energy is made in both pure natural single crystals 
of graphite and in several highly ordered pressure annealed samples of pyrolytic graphite. 
This is done by observing splittings of Landau level crossings of the Fermi energy as the 
magnetic field is increased. Fields ranged to 10 teslas and temperatures were near 
1. 1 K. Quantities studied were the electrical magnetoresistance, the Hall effect, and the 
thermoelectric power. To interpret the data in terms of existing theories, shifts in the 
Fermi energy with magnetic field were studied. Experimental splittings are compared 
with a theory by McClure and Yafet and discrepancy is partially accounted for by large 
shifts in Fermi energy with magnetic field. 


INTRODUCTION 

Graphite is an important aerospace material because of its very high Young’s modu- 
lus, and its extremely anisotropic and unusual properties (ref. 1). It has been formed 
into composite material for high strength applications, for example, reference 2. Graph- 
ite is a layered structure of carbon where the in-plane covalent bonding strength is as 
great as that for diamond. Highly perfected graphite (pyrolytic graphite, ref. 3) has only 
recently been synthesized, so it is not well studied or understood. The present report 
discusses one of many discoveries made in graphite while studying the very high mag- 
netic field properties. These properties were found to be extremely anisotropic and thus 
potentially very useful in devices. More high field graphite work is reported in refer- 
ences 4 to 9. 

The present report presents the first discovery and study of the effects of electron 
spin on the electronic structure of graphite in a strong magnetic field. It presents new 
evidence that pyrolytic graphite and perfect natural single crystals have much in common. 



This report shows that the parameters needed to characterize other properties of pyro- 
lytic and crystalline graphite will have to be modified to account for the spin effects re- 
ported herein. 

In the presence of a very strong magnetic field, a conducting solid at low tempera- 
tures will have an electron distribution in momentum space characterized by Landau 
levels. Landau levels are cylinders with axes parallel to the magnetic field and contained 
within the zero field Fermi surface. Each cylinder corresponds to a discrete energy 
level for the electrons. Each Landau level is described by an orbital quantum number n 
and a spin quantum number having two values: one for spin up and one for spin down. In 
most materials the spin energy levels are degenerate. In graphite these levels are very 
slightly separated. There is a substantial spin angular momentum - orbital angular mo- 
mentum coupling which causes ’’splitting" of the normally degenerate spin levels. The 
subject of this report is the study of the strength of the spin-orbit coupling. To do this 
properly requires a study of the energy shifts as a function of magnetic field. Fortu- 
nately both the spin splitting and Fermi energy shifts can be obtained from the same ex- 
perimental data. This assures that the sample and experimental conditions are identical 
for both studies. Spin splittings are observed as double extrema in several transport co- 
efficients as a function of field, and these are a measure of the conduction electron g- 
factors when shifts in the Fermi energy with field are accounted for. 

Wagoner (ref. 10) observed conduction electron g-shifts from electron spin resonance 
studies. Spin resonance, however, measures the g-shifts averaged over the Fermi sur- 
face for particular magnetic field directions. Splittings of extrema in transport oscilla- 
tions due to spin split Landau levels measure the spin-orbit splitting for particular mag- 
netic field directions. Splittings of extrema in transport oscillations due to spin split 
Landau levels measure the spin-orbit splitting for particular points in the Brillouin zone 
(refs. 11 and 12). Knowing the location results from being able to identify a series of 
Landau level crossings of the Fermi energy Ep with a known location of the particular 
Fermi surface in the Brillouin zone (refs. 4 to 6, 9, and 13). 

Figure 1 shows the Brillouin zone and Fermi surface of graphite. The Brillouin zone 
is hexagonal with the C-axis perpendicular to the hexagonal plane, as shown in figure 1. 
The Fermi surfaces are located along the vertical edges of the zone. Sections A are hole 
carriers and sections B are electrons (refs. 4 to 6, 9, and 13). McClure and Yafet 
(ref. 14) have calculated the g-shift for particular locations in the Brillouin zone. How- 
ever, to compare with the spin resonance experiment, they had to take an average over 
the electron and hole Fermi surfaces in graphite for the magnetic field along the C-axis. 
This report presents results which can be compared directly with the McClure and Yafet 
theory, which in turn is based on the Slonczewski -Weiss (SW) band model (ref. 15). 
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EXPERIMENTAL RESULTS 


The magnetoresistance is measured by a conventional four-probe technique and tem- 
peratures were regulated near 1. 1 K by pumping on a liquid helium bath. For more de- 
tails of this method see references 16 and 17. 

Figure 2 shows the magnetoresistance of a natural single crystal as a function of 
magnetic field to 10 teslas (1 tesla = 10 kG) at 1. 1 K for the magnetic field parallel to the 
[0001] axis (C-axis). Adams and Holstein (ref. 18) have shown that conductivity maxima 
occur at coincidence of Landau levels and Ep. We experimentally find (refs. 4 to 6, 
and 9) that in graphite this implies resistivity minima at Landau level crossings of Ep. 
References 6 and 9 show why this is true. Essentially it resvilts from the Hall resistivity 
being much smaller than the magneto resistance. 

In figure 2 there are distinct minima at 6. 60+0. 02 and 7. 40±0. 02 teslas, where the 
n = 1 Landau level for electrons (in refs. 4, 5, and 13 it is shown that the electrons are 
located at B in fig. 1) crosses the Fermi energy. In highly ordered pressure annealed 
samples, we find resistivity minima less pronounced than for single crystals but occur- 
ring at nearly the same fields as for single crystals. An example of the pyrolytic data is 
shown in figure 3, where Hall resistivity p,,„ and Hall conductivity a multiplied by 

2 ^y 

H are plotted, as well as the resistivity as a function of magnetic field at 1. 1 K. 
The thermoelectric power as a function of magnetic field is shoAvn in figure 4 for a tem- 
perature of 1. 1 K. In pyrolytic samples we observe spin split maxima in the Hall con- 
ductivity (fig. 3) and minima in the thermoelectric power for the n = 1 electron level. 

The field values for extrema and splittings in resistance. Hall effect, and thermopower in 
natural single crystals as well as in several pyrolytic samples are listed in table I. In 
single crystals, figure 2 shows the spin split n = 1 hole Landau level resistivity minima 
appearing at 3. 643+0. 01 and 3. 412+0. 007 teslas. The spin split n=2 electron level 
crossings appear at 2. 936+0. 005 and 2. 992+0. 007 teslas. 


THEORY 


The energy bands in zero field along the vertical edge of the hexagonal Brillouin zone 

edge are shown in figure 5. (See ref. 19 for a brief description. ) The distance, in the 

k k plane, from the vertical edge is a. For nonzero the energy eigenvalues are 
X y 


E = 1 (Ej + Eg) ± |1 (Ej - Eg)2 + y2(i - ,.) (1) 

E = A (E 1 + Eg) + I^A (Eg - Eg)2 + yg( 1 + (2) 
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which come from solutions of the effective mass Hamiltonian (refs. 15 and 19). (All 
symbols are defined in the appendix. ) In equations (1) and (2) 


Ej= A + y^r + -r5r 


1 - -2 
2 


E2= A-rir + ^rgr^ 


^3 = ^ 2 ^^ 


(3) 

(4) 

(5) 


and 


r = 2 cos 7T| 


( 6 ) 


where ^ is the normalized distance along the vertical edge of the Brillouin zone shown 
in figure 1. In figure 5, a = 0 and equations (1) and (2) reduce to Ej, E2, and Eg. In 
the preceding equations 


The band parameters are 



(7) 


Yq = 3. 13 [2. 85] eV 
yj = 0.4 [0. 31]eV 
Y 2 = -0.0183 [-0.0185] eV 


yg = +0. 3 [0. 29] eV 
y^ = 0. 25 [0. 18] eV 
yg = -0. 0183 eV 
A= 0,025 eV 
a^ = 2. 456X10"® cm 




( 8 ) 
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which are determined by experiment. (See ref. 19, for example). The bracketed values 
in equation (8) and in equations which follow are from magnetoreflection data, unchecked 
against de Haas -van Alphen results (ref. 20). 

The energy levels in a magnetic field were worked out by McClure (ref. 21) and Inoue 
(ref. 22), and the resultant secular equation is 



E - Eg 

E “ ^ ” ^2 

4« 4 

fE -Esp 

E - Ej 

_ E - Ejl 2 ^ q2'_ 

2 



L ^ L( 


(1 + pA ^ 


(9) 


and 


Q = 


3a2y2 eH 

o^O 


2Rc 


( 10 ) 


Constants e, R, c have their usual meaning, and H is the magnetic field. Equation (9) 
thus relates the energy E to the magnetic field strength H for particular v values de- 
termined by equations (6) and (7). 

The energy levels in a magnetic field can be found approximately by (e.g. , ref. 11) 



( 11 ) 


( 12 ) 


(13) 


This equation is used experimentally to obtain the g-factor. 

McClure and Yafet (ref. 14) have worked out g-factors and shifts 6g from the free 
electron value of 2. 002, due to spin-orbit coupling. They get (after correction for a 
minus sign) 
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( 14 ) 



2 2 r 2 2 1 

3Aa yQm^|2A + Si/y^cos + 4(yg - ^ 2)^08 7 r|J 
4R^ I Eg - EpjD^yj^cos ?r| 


where 


^ _ 2Ep - (Eg + Eg) 

+ 2 
(1 + i; cos 7T I) 

D = ~ ^^1 + ^3^ 

2 

(1 - y cos JT^) 


(15) 


(16) 


where v is given in equation (7) and the plus signs are for electron carriers. Param- 
eter X is the spin-orbit coupling constant and is a measure of the strength of the spin- 
orbit interaction. Other parameters and constants are the same as discussed earlier. 
Parameter yg does not appear in any of the preceding equations because it was neglected 
in these derivations to simplify the analysis. Equation ( 8 ) shows that the magnitude of 
yg is large compared to several other constants. Since these parameters enter the ele- 
ments of the Hamiltonian linearly (ref. 23), this is a serious neglect. 


Interpretation and Analysis of Experimental Results 

In figure 2 the resistance minima determine the magnetic field values for crossings 
of Ep and the Landau levels. (This is discussed in detail in refs. 4 to 6 , and 9). Using 
the experimental minima from figure 2, equation (9) is solved for E by computer. The 
resultant E values are then adjusted (by subtracting 2yg) such that they become Fermi 
energy values. These numbers for electron carriers are 


E = +0.0230 [0.0231] eV 

n= 1 

(17) 

E = +0.0169 [0.0170]eV 

n = 2 

(18) 

E = +0. 0156 [0. 0158] eV 

n = 3 

(19) 

(The bracketed values are from ref. 19. ) Using a similar analysis 

for the hole carriers 


gives 
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E = 0.0164 eV n = 1 


( 20 ) 


In figure 6 the energy values from equations (17) to (20) are plotted as a function of 
the magnetic field associated with the n = 1, 2, 3 Landau level crossings. The mean 
value of magnetic field is taken for the n = 1 levels where the splitting is large. For the 
n=l electron spin split Landau level crossing, the values of and E_ are found 
from the graph. The field values for crossing the Fermi energy for spin+ and spin- 
are known from figure 2 (for the natural single crystal) and listed in table I. The curve 
in figure 5 thus gives 


E^(6. 60 T) = 0. 0222±0. 0002 eV (21) 

E_(7. 40 T) = 0. 0237+0. 0002 eV (22) 

for the n = 1 electron Landau level. From equation (13) this gives 

g = 2. 7+0. 4 (23) 


so 


5g = 0. 7±0. 4 n = 1 electron (24) 

For the n=2 electron level E^ and E_ are approximately 0. 0170±0. 0003 eV and 
equation (13) gives a g such that 

6g = -O.OiO.3 n = 2 electron (25) 

For the n = 1 hole level, E is 0.0164 and equation (13) gives a g such that 

6g = 3. 0+0. 5 n = 1 hole (26) 

However, if the Fermi level shift is accounted for as shown in figure 5, then 

E^=0. 0164eV (27) 


and 


E = 0.0168 eV 


(28) 
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This means that 


6g = +1.4 (29) 

McClure and Yafet analyzed Wagoner’s data (ref. 10) using equation (14) and obtained 

A ~ 3x10"^ eV (30) 

However, was assumed to be zero, Permi energy shifts with magnetic field were not 
considered, and a possible numerical error exists (J. W. McClure: private communica- 
tion). Dresselhaus and Dresselhaus (ref. 24) have calculated the spin-orbit interaction 

in graphite, avoiding some of the simplifying assumptions made by McClure and Yafet, 

-2 

and find that X should be on the order of 10 eV. Using McClure and Yafet’s theory and 
the experimental g- values, we get 


X = -1.4+0.8x10"^ eV 

n = 1 electron 

(31) 

A = 1. 4+4x10"^ eV 

n = 2 electron 

(32) 

A = 2.3+0. 2x10"^ eV 

n = 1 hole 

(33) 

These equations indicate that the magnitude of A 

-3 

is on the order of 10 eV, 

which is an 


order of magnitude greater than McClure and Yafet’s value but is smaller than found in 
reference 24. The values found in equations (31) to (33) do not agree with each other. 

The solution to this problem is discussed in the next section. 

DISCUSSION AND CONCLUSIONS 

Thus, in this report, a study of the spin-orbit coupling energy is made in both pure 
natural single crystals of graphite and in several highly ordered pressure annealed sam- 
ples of pyrolytic graphite. This is done by observing splittings of Landau level crossings 
of the Fermi energy as the magnetic field is increased. Fields ranged to 10 teslas and 
temperatures were near 1. 1 K. Quantities studied were the electrical magnetoresistance, 
the Hall effect, and the thermoelectric power. To interpret the data in terms of existing 
theories, shifts in the Fermi energy with magnetic field were studied. 

In order to find the experimental g-shift using equation (13), the Fermi energy at 
spin-up and spin-down Landau level crossings has to be found from equation (9). Equa- 
tion (9) is based on the SW band model and depends on the parameters of equation (8), 
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which need to be better determined. Equation (9) also neglects y^, which may signifi- 
cantly affect the numbers found in equations (17) to (20) and thus the experimental g- 
value. This is probably the cause of the sign discrepancy with the n = 1 electron g- 
value. 

Equation (14) for the theoretical g-shift also neglects y^, makes other simplifying 
assumptions, and depends on accurate values for the constants of equation (8). Thus 
there is a need for a determination of the best band parameter constants consistent with 
all available graphite data. There is also a need for more theoretical work to predict the 
relation between g-shifts and spin-orbit coupling energy in graphite, in order to better 
compare experiment with theory. 

McClure (private communication) has recently pursued this theoretical work, in- 
cluding effects of yg and Fermi level shifts, and gets approximate agreement between 
theory and experiment. The results depend critically on the parameters of equation (8). 

J. Furdyna (private communication) has observed the effects of spin splitting of the 
n = 1 electron level on helicon propagation in a sample of pyrolsdiic graphite. The split- 
ting appears roughly at the fields discussed in this report, for the n = 1 electron level. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 4, 1970, 

129-02. 
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APPENDIX - SYMBOLS 


a,a^ 

B 

c 

D^.,D_ 

E 

E„ 


± 

e 


g 

ag 

H 

H+,H_ 

h 

m* 


m, n 

“o 

Q 

Of 


y 

>b’ 

y ^, 75, ^ 

r 

V 

a 

X 


locations of holes in Brillouin zone 
graphite lattice parameter 
locations of electrons in Brillouin zone 
velocity of light 

parameters in McClure, Yafet theory 
energy 

Fermi energy 

energy relative to band edge for spin up and spin down 

parameters in Slonczewski -Weiss (SW) band model 

electric charge 

spin splitting factor 

spin splitting factor shift from 2. 002 

magnetic field strength 

field values for spin up and spin down 

Plancks constant divided by 2 tt 

effective mass 

integers 

mass of free electron 

parameter in SW band model in a magnetic field 
field dependent parameter 
phase factor 

constants of SW band model 

parameter in SW band model 

parameter in SW band model 

dimension along vertical edge of Brillouin zone 

distance from vertical edge of Brillouin zone in plane 

spin-orbit splitting constant 
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TABLE I. - FIELD VALUES FOR SPIN SPLIT LANDAU LEVEI£ 


Sample 

Landau level 

Magnetic field values, 



Magneto- 

resistance 

Hall 

coefficient 

Thermopower 

Natural 

single 

crystal 

n = 1 electron 

7. 40±0, 02 
6. 60±0.02 

2. 992±0.007 
2.936±0.005 





n = 2 electron 






n = 1 hole 

3. 643±0.01 
3.412±0.007 








Pyrolytic 
graphite:^ 
PG 5 

n = 1 electron 

7. 54±0.05 
6.92±0. 1 

7. 56±0.05 
6. 93±0. 1 


PG 4 

n = 1 electron 

7. 54±0.05 
6. 92±0. 1 

7. 53±0.05 
6.82±0. 1 

7. 68±0.05 
6.92±0. 1 

PG 3 

n = 1 electron 

7. 23±0. 2 
6. 43±0.2 

7.45±0. 1 
6.80±0. 1 



T = 10 kG. 

*^PG 5, PG 4, and PG 3 refer to different samples of highly 
ordered pressure annealed pyrolytic graphite. 


(C) 


R 




Minority carrier 
pockets predicted 
by SW theory^ 



Figure L - Brillouin zone and Fermi surfaces of graphite. C-axis direction is 
marked. Section A are holes; section B are electrons; H and K are points 
along the Brillouin zone edge. (From ref. 19. ) 


Magneto resistance, A p/p 



t. 


60 


n = 1 electron 


Figure 2. - Magneto resistance as function of magnetic field at 1. 1 K. Field values at spin split minima are in teslas 
(1 T = 10 kG). Natural single crystal graphite. 


n = 3 

electron^ 


n * 1 electron: 



Figure 3. - Magneto resistance pyy, Hall resistivity, and Hall conductivity 
times H^, as function of magnetic field at 1. 1 K. Spin splitting of n = 1 
electron levels are marked. 





igure 6. - Fermi energy as function of magnetic field, as obtained from experimental data 
of figure 2 and equation (9). 
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